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Abstract 

We study medium modifications of the the e + e~ and/or p + p~ mass spectrum in 
coherent photo- and electroproduction of the p°(lS)- and p'(2S')-meson resonances on 
nuclear targets. The analysis is performed within the coupled p°(lS), p'(2S), ... channel 
formalism in which nuclear modifications derive from the off-diagonal rescatterings. We 
find that the effect of off-diagonal rescatterings on the shape of the e + e~ mass spectrum in 
the p°( 15*) -meson mass region is only marginal but is very important in the p'{2S) mass 
region. The main off-diagonal contribution in the p'(2S) mass region comes from the 
sequential mechanism 7* — > p°(lS) — > p'(2S), which dominates the p'{2S) production for 
heavy nuclei. Our results show also that in the p'(2S) mass region there is a considerable 
effect of the interference of the Breit-Wigner tail of p°(lS')-meson with the p'(2S')-meson. 



1 Introduction 



The in-medium modification of hadrons in the cold nuclear matter and hot QCD medium 
has been under active investigation during the last years. In particular, there was con- 
siderable interest in the medium effects for light vector mesons. In several papers the 
masses of vector mesons at rest in nuclear matter have been calculated within different 
approaches [M 0, BL 01, EjL 0, M . The in-medium effects for moving vector mesons have been 
discussed in refs. || ||, |10|, [ll], HH- Recall that in optics medium effects are described by 



the refraction index which for the dilute media is calculable in terms of the photon-atom 
forward scattering amplitude. Extension of this formalism to fast vector mesons with the 
wavelength A = 1/p much shorter than separation of nucleons in the nuclear matter gives 
the in-medium mass shift and collisional broadening the form (||, for the related early 
works see 113, [bl) 



Am v (E) = -27i^Ref(E) , AT V {E) = —pa{E) , (1) 
my THy 

where E and my are the meson energy and in- vacuum mass, f{E) is vector meson- 
nucleon forward scattering amplitude and ha is the nuclear matter density. These formulas 
are quite general and must hold for any particle in an infinite medium if the inelastic 
rescatterings and/or coupled-channel effects can be neglected, see below. An experimental 
observation of the in-medium mass shift and collisional broadening (1) would be very 
interesting. Potentially it could give information on the VN scattering amplitude which 
cannot be measured directly. 

In order for the mass shift and collisional broadening (1) to be observed experimentally 
the typical decay length L d ~ E/m v Ty must be smaller than the nucleus radius Ra, i.e., 
the momentum of the vector meson must be smaller than 

p p < 6 GeV, p w < 300 MeV, p^ < 200 MeV , (2) 

for the p -, uj -, and 0°-meson, respectively ||. Thus only for the p°-meson [] there is a 
sufficiently broad energy interval where the relations (1) could be used. 

The applicability limits (§) are purely kinematical ones and do not take into account 
possible interference of decays of the vector meson inside and outside the target nucleus 
and the quantum effects in production of the vector mesons. The both effects were 
studied, and found to be important, in the recent |H| Glauber-Gribov multiple-scattering 



theory [|15|, [16! analysis of coherent p°-meson photoproduction observed in e + e~ and/or 
p + /!~ dilepton mode. The dilepton mass spectrum was shown to have a two-component 
structure corresponding to the decays of p°-meson inside and outside the target nucleus. 
The inside component can be approximately described by a Breit-Wigner formula with 
the in-medium modified mass and width as predicted by Eqs. (1). However, since the 
nucleus has a finite size, this component with a broad width does not develop a genuine 
pole in the complex plane of the invariant mass, M, of the lepton pair. The genuine 
pole in the complex M-plane pole comes only from the p°-meson decays in vacuum and 
this outside component can be described by the standard Breit-Wigner formula with 

1 Hereafter wherever appropriate the p° will stand for the ground state p°(lS) meson. 
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the in-vacuum mass and width. Ref. \W\ found that even at low energy E ~ 2 GeV 
such that Ra/L^ ~ 2 for heavy nuclei, the interference between the inside and outside 
components is substantial and produces a complex dilepton mass spectrum which cannot 
be described in terms of a Breit-Wigner formula with a definite mass and width. For 
instance, at E — 2 GeV the dilepton mass spectrum was found to develop a minimum 
near the p°-meson mass. An experimental observation of this phenomenon would be of 
great interest. On the theoretical side, this calls for a numerical analysis within tested 
models of photoproduction on how a model analysis of such data would allow to distinguish 
the inside and outside components and extract the p°-meson mass shift and collisional 
broadening for the inside component. 

In the present paper we study the coherent reaction j*A — > VA — > e + e~A,p + p~A, 
where 7* is a real or virtual photon and the target nucleus remains in the ground state, 



within the coupled-channel approach extending over the early work |T(J in which the 
numerical calculations have been performed neglecting the off-diagonal rescatterings of 
p°-meson. The coupled-channel analysis presented here is based on our well tested colour 
dipole approach (see, for instance, |T7| and references therein), which was earlier success- 
fully used in the analysis of the data on p° and J/\l/ electroproduction on nuclear targets 



at high energies [H], |2(| and vector meson production at HERA [fL7|, f21| , This 



analysis is of interest for two reasons. First, including the off-diagonal rescatterings allows 
one to check the accuracy of the one-channel approximation in the p°(lS')-meson mass 
region 0.5 < M < 1 GeV studied in ref. |TIJ. On the other hand, in the coupled-channel 
approach one can extend the mass region and investigate the medium effects for the 25' 
state p'{2S) for which the sequential off-diagonal mechanism 7* — > p°(lS) — > p'{2S) is 
potentially important. This extension to the p'(2S')-meson is of great interest in itself. 
The key feature of photoproduction of the 25" vector mesons on a free nucleon is strong 
suppression due to the nodal structure of the wave function of the 2S state [f23], [T8J. In ref. 
fl2~4j it was shown that the node effect can lead to an anomalous A- and Q 2 -dependence 
of p'(2S) photo- and electroproduction. This effect may help to resolve the long standing 
problem of the D-wave vs. 25-wave assignment for the p'(1480) and p'(1700) states. The 
strong suppression of the cross section for the p'(25')-meson as compared to the p°(lS')- 
meson makes the experimental study of this phenomenon a challenging task. In particular, 
the mass spectra of the final particles in the p'{2S) mass region can be affected by the 
interference with the p°(lS')-meson Breit-Wigner tail. Our colour dipole coupled-channel 
approach describes very well the suppression of ty' production as compared with J/\l/ 
observed by the NMC [B^, 26], E687 collaborations and HI collaboration pS|| and 
provides a sound framework for understanding the the prospect of experimental study 
of the p'(2S')-meson production 7*v4 — > p°(lS)A, p'{2S)A — > e + e~A. In the reported 
numerical analysis we focus on the range of energies of the forthcoming high-luminosity 
experiments at TJNAF . 

The paper is organized as follows. In section 2 we give the basic formulas of the 
coupled channel approach to the coherent j*A — > e + e~A,p + p~A reaction. In section 
3 we discuss evaluation of the diffraction scattering matrix. The numerical results on 
nuclear modifications of the in-nucleus decay component of the dilepton mass spectrum 
are presented in section 4.. The results are summarized in section 5. 
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2 The coupled channel formalism 



The coupled-channel formalism for the 7 A — > e + e A reaction is found in ref. (10 



sec 



also our early work |ZIJ . For this reason, here we discuss it briefly and give only the basic 
formulas which are necessary for understanding technical aspects of our approach. 

The standard Glauber-Gribov multiple-scattering theory [|15|, [16|] for coherent inter- 
action of the projectile of energy E with a nucleus is equivalent to solving the set of 
coupled-channel eikonal wave equations [ fj| : 







•2 



dz 2 



+ m 2 + U(r) 



^•(r) = £ 2 ^(r). (3) 



Here the z-axes is chosen along the photon momentum; is the wave function for the 
channel |z) which can be a hadronic resonance or the initial photon; m 2 is the diagonal 
mass operator with eigenvalues m 2 = (mf — imiT^Sij, where rrii and Tj are the in- vacuum 
mass and width of the state for the incident photon T 7 . = and m 2 * = — Q 2 where 
Q 2 is the photon virtuality and the optical potential U(r) in (^) equals 

U l3 {v) = -A7t(i\f\ J )n A {T) ) (4) 

where n^(r) is the nuclear number density and / is the forward scattering matrix in the 
normalization 

Im(i|/|i) = ^a tot {iN -+ iN) . 
Air 

The boundary condition on the front face of the nucleus, z = —Ra, is \i) = 8^* . 

The probability amplitude for the coherent transition ^*A — > e + e~A be written in 
terms of the solution ^(r) of eq. (|]) as fpTOfl 

T(E,M,p ± ) = NY,(e + e'\t\i) J d 2 bdzexp[-i(p z z + p ± b)]^(r) , (5) 

i=h 

where r = (b, z), M is the invariant mass of the e + e~ pair, p z and p ± are its longitudinal 
and transverse momenta, respectively, (e + e~\t\i) is the probability amplitude for % — > e + e~ 
transition, is a normalization factor which is immaterial from the point of view of the 
shape of the e + e~ mass spectrum. The summation in (||) goes only over the hadronic 
states, which in our case are vector mesons p°(lS), p'(2S) etc. For a heavy nucleus with 
a mass number A ^> 1 and at E ^> M the nuclear recoil can be neglected and the 
longitudinal momentum of the e + e~ pair equals 

Q 2 + M 2 + pj 
Pz 2E ' ( ' 

In the coherent production < 1/R\ and in what follows we focus on p ± = and 
suppress this argument. Then by virtue of @ the p 2 -dependence of the spatial integral in 
the right hand side of Eq. (g) transforms directly into the M-dependence of the amplitude 
T(E, M, p_|_) and nuclear modification of the e + e~ spectrum. 
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At high energies the solution of (^) for the hadronic sector which we need for evaluation 
of the amplitude (^) can be written in the form 



V h (b,z) = (h\S(b,z)\^*}exp(ip r z), 
where the operator S is given by 



S(b, z) = P z exp 



2p, 



d£ m 2 + Q 2 + U(b,g) 



(7) 



(8) 



Here, P z is ^-ordering operator, and p y * is the photon momentum. For the numerical 
calculations it is convenient to treat in Eq. ([8]) the off-diagonal part of the optical potential 
in the hadronic sector as a perturbation while the diagonal transitions are included to all 
orders. Then, following Ref. one can represent the matrix elements of the operator 
S(b, z) in the form of the i^-fold off-diagonal rescatterings series 



(/»|^(b^)|7*> = E<^ M 0>^)l7*> 



u=0 



where 



(9) 



I f 1 

(h\S {Q) (b,z)\~j*) = --a hr J dz 1 n A (b,z 1 )exp[ik J * h (z-z 1 )--t(b,z,z 1 )cr hh \, (10) 



(h\S {u) {b,z)\^*) = (-- J £ a hii/ a' ii/ivl ■ ■ ■ a il7 * exp(ik 7 * h z) 

f r 1 

x J dz u+1 n A (b, z u+l ) exp[ik hiv z u+1 - -i(b, z, z u+x )a hh 

— oo 

x J dz v n A (b, z v ) explik^i^Zv - -t(b, z u+1 , z v )a iviv \ ■ 

— oo 

22 1 
x J dz 1 n A (b, z x ) exp[zfc il7 .zi - -t(b, z 2 , zi)a hh ] , 



v > 1. 



(11) 



Here, a ik = — Sik&M, the matrix a is connected with the forward diffraction scattering 
matrix 



Z2 



t(b,z 2 ,z l )= I dzn A (b,z) 
is the partial optical thickness, and 



k; 



kh-y* — k. 



mf — irriiTi — m 2 + irrijTj 

2E : 
m 2 - irriiTi + Q 2 



7*/i 



IE 



(12) 
(13) 

(14) 
(15) 



The exponential factor (z — z\) — jt(b, z, zi)au] in (|i~0D, 

onal, iN rescatterings to all orders. 

The real part of (|T5|) - the longitudinal momentum transfer in j ' i transition - rises 
with the difference between rrii and rrij. Consequently, the oscillating exponential factors 
in (|I0D, (p])) lead to the form factor suppression of the contribution from heavy interme- 
diate resonance states and the related form factor suppression of heavy mass production 
and of the coherent cross section at large Q 2 , when the longitudinal momentum transfer 
k^* in the 7* —>■ h transition becomes large. Precisely the same form factor effect gener- 
ates strong dependence of the probability amplitude (H) on the mass of the e + e~ pair so 
that the shape of resonances would differ strongly from the standard Breit-Wigner one. 

Eq. (H) quantifies the separation of the production amplitude T(E, M) into the inside 
and outside components. Namely, beyond the target nucleus the z-dependence of wave 
functions ^i(z, b) follows the in- vacuum decay law and the corresponding outside con- 
tribution to T(E, M) has the familiar form of the sum of Breit-Wigner amplitudes with 
residues proportional to ^i(z = -hR^b). The nuclear effects cam modify dramatically 
both the relative amplitude and phase of these residues compared to a hl * for the free 
nucleon case. As explained in the Introduction, the inside contribution is a Fourier trans- 
form over the finite range —Ra < z < Ra and could develop the Breit-Wigner form with 
the shifted mass and collision-broadened width only provided that La <C Ra- 

The first order term (|i~0|) corresponds to the standard Glauber approximation when the 
state \h) produced in the 7* — > h transition then propagates through the nucleus without 
inelastic rescatterings. The correction from the off-diagonal rescatterings is given by (]TTf) 
and is responsible for the color transparency phenomenon in electroproduction of vector 
mesons at high Q 2 where it changes drastically the cross section as compared with the 
Glauber model predictions. However, as will be seen from our results for — > p'(2S)A 
reaction the off-diagonal effects come into play already at moderate values of the photon 
virtuality Q 2 < 1 GeV 2 . This is a consequence of strong suppression of the direct 7* — > 
p'{2S) transition. As a result, the sequential mechanism 7* — > p°(lS) — > p'(2S), involving 
the off-diagonal p°(lS) — > p'(2S) rescattering, becomes important even at Q 2 = 0. 



IT) sums elastic, diag- 



3 Calculation of the forward diffraction matrix 

As an input to the coupled- channel calculations one needs the forward diffraction matrix. 
At high energies it can be written as a sum of the Pomeron and reggeon contributions 

a = a-jp + a R . (16) 

At energies E > 2 GeV to be considered in our paper, the dominating contribution to 
the a comes from the Pomeron exchange. We evaluate the Pomeron exchange within the 
dipole approach describing the resonances as nonrelativistic qq states. Then the Pomeron 
contribution to the diffraction matrix element cjjj for hadronic states can be written as 

(i\ajp\k) = J d 2 pdz^/j*(p, z)a(p)ip k (p, z) , (17) 

where p is the transverse size of the qq pair, ipi^p, z) are the wave functions describing 
the qq states, and a(p) is the cross section of interaction of the qq pair with a nucleon. 
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We also need the excitation matrix elements (z|<t|7*) for the 7* — > qq excitation on 
a free nucleon which in the non-relativistic approximation proceed into qq states with 



the sum of quark helicities equal to the photon helicity [23]. Using the corresponding 



perturbative light-cone wave function of the virtual photon ||30|| , we have |2 

(z|<7| 7 *> = C J d 2 p^( Pl z = 0)a(p)K o (ep) , (18) 

where 

e 2 =mJ + Q 2 /4, (19) 

m q is the quark mass, K (x) is the modified Bessel function. In the present paper we 
focus on the nuclear-modification of the shape of resonances, and the absolute value of 
the normalization factor C in Eq. ([18]), and in Eq. ([|) as well, is immaterial. 

We use the oscillator wave functions for the qq states, which simplifies considerably 
the numerical calculations: because of azimuthal symmetry of crjp(p) only the off-diagonal 
rescatterings change only the transverse excitations with zero azimuthal quantum number 
can be excited in the intermediate state. Excitation energy of the transverse qq oscillator 
equals 2hu, here u is the oscillator frequency. We take in our analysis the quark mass 
m q = rn p /2 and the oscillator frequency uj = [ra p i — m p o)/2 « 0.35 GeV , assuming 1480 
MeV for the mass of the radial excitation of p°-meson. For the widths of the first two 
excitations we use the values r p o w 150 MeV and F p / w 285 MeV and, following the string 
model |3l], assume Tj cx m, for higher lying states. We are not sensitive to this latter 
assumption, though, because the contribution of the higher excitations turns out to be 
very small. 

The experimental data on the low- 2 structure function F 2 and vector meson elec- 
troproduction off nucleon can be described by representing a(p) as a sum of the energy 
dependent perturbative and energy independent nonperturbative components |i~7|, |22|, |32|| . 
Here we focus on the energy region E < 20 GeV where the energy dependence of the dipole 
cross section can be neglected to first approximation. Also, in the present analysis we are 
sensitive to a(p) mostly in the nonperturbative region of p > 0.5. Here the gross features 
of a(p) are well parameterized by the two-gluon exchange model of the Pomeron |33}| : 



a{p) = 1M / rf> q [l-exp( i qp)]|l-G 2 ( q ,- q )] 
3 J (q 2 + p 2 g ) 2 

here G2(qi, q2) = (^1 exp(iq 1 ri + zq 2 r 2 )|iV) is the two-quark form factor of the nucleon, 
Pg0.3 GeV is an infrared cutoff. It reproduces the colour transparency property <r(p) oc p 2 
at small p, by which the point-like qq system, which can be represented as a superposition 
of an infinite set of the resonance states, propagates through the nucleus without interac- 
tion. This color transparency phenomenon in terms of the resonance states is connected 
with exact cancellation of the diagonal and off-diagonal amplitudes []35|| . The coupling 
constant as was normalized so the Pomeron contribution to the p°N total cross section 
<jp(p°N) = (p |o"ip|p ) ~ 20 mb. In this case the two-gluon formula gives at p > 0.5 
fm the a(p) which is close to the dipole cross section extracted from the analysis of the 



experimental data on electroproduction of vector mesons |T7j . Hereafter we shall consider 
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real, Q 2 = 0, and virtual, 
functions we find 

R(2S/1S) 



Q 2 = 1 GeV 2 , photoproduction. With our nonrelativistic wave 



(p'(2S)\a\r) 




Q 2 
Q 2 



o. 

1 GeV 2 



(21) 



(p°(lS)\a\r) 

which is close to predictions of ref. obtained for the relativized wave functions. At 
p < 0.5 fm the present parameterization gives a(p) somewhat larger than that of ref. 



17| . However, because of the larger quark mass in the present nonrelativistic model for 
qq states the resulting diffractive matrix turn out to be close to that of ref. [|I7|]. At 
this point, it must be made clear that the two-gluon parameterization (|20"D of cr(p) is 
oriented towards description of the combined nonperturbative+perturbative dipole cross 
section and p g is a phenomenological parameter which must not be taken at the face 
value. The analysis of low-x HERA data on the proton structure function F 2p within 
the generalized BFKL equation [36], and the nonperturbative evaluation of the gluon 
correlation radius [37] yield a clearcut evidence in favor of the infrared cutoff pc ~ 0.75 
GeV for the perturbative dipole cross section. 

In our calculations we take into account the first 4 transverse excitations. The Pomeron 
contribution to the diffraction matrix in terms of these transverse oscillator states which 
we obtained from Eq. (|i~7"|) for our parameterization ( ^0|) of a(p) is given by 

( 



(i\ajp\k) 



Here the matrix elements in 



V 



20 


-10.1 


-4.8 


2.4 \ 


-10.1 


30.2 


-8.0 


5.2 


-4.8 


-8.0 


32.9 


6.5 


2.4 


5.2 


6.5 


34.1 ) 



(22) 



ik 



|) are in units of mb, i and k are the radial quantum 
numbers of the transverse oscillators (as was said above the Pomeron exchange does not 
change the longitudinal quantum number). The matrix (|22|) shows clearly the decrease 
of the off-diagonal amplitudes with increase of the difference between the initial and fi- 
nal radial quantum numbers \i — k\ which derives from the oscillation of the resonance 
wave function and, in conjunction with the form factor effect, suppresses the contribu- 
tion of higher excitations in the production amplitude matrix element (O). Using the 
three-dimensional p'(2S")-meson wave function from ( p2|) one obtains for the Pomeron 
contribution to the p'(2S)N total cross section a$(j/(2S)N) = (p'(2S)\ajp\p'(2S)) « 27 
mb. 

In parameterizing the reggeon contribution to the diffraction matrix we assume that 
the secondary reggeon exchanges can be treated in terms of scattering amplitudes for the 
quark (or antiquark) making up the qq state as predicted by the dual parton model P| 
and quark gluon string model |39| based on the idea of the topological expansion In 
this case one can neglect the reggeon contribution to the off-diagonal transitions, and for 
all excited p' states the reggeon contribution to diagonal p'N scattering amplitudes turn 
out to be equal to the reggeon contribution to the p°(lS)N scattering amplitude. This 
amplitude is dominated by the contribution of the Regge pole P' which can be written in 
the Regge approach as 



(pVrIp ) = pp> 



so 



a pi — I 



1 + cos7rap/ 
sin7rap/ 



(23) 
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where s = m 2 N + 2ErriN. In our analysis we take the standard reggeon intercept api = 0.5. 
The residue rp> has been adjusted to reproduce at E ~ 10 GeV the real part of the 
p°N scattering amplitude extracted in ref. ||T^] from the experimental data for p°-meson 



photoproduction using the vector dominance model. For sq = 1 GeV 2 this gives rpr « 15 
mb. The nonzero Re/Im ratio for the diffraction scattering matrix leads to the mass 
shift for the resonance states decaying inside the nucleus. For the p°-meson in the energy 
region E ~ 2 — 20 GeV considered in the present paper it gives Am p o ~ 50-100 MeV. 

For evaluating the e + e~ mass spectrum we need also the transition amplitude 
(e + e~\t\i) which enters Eq. (||). Following the authors of ref. [|Kj we neglect a pos- 
sible smooth M-dependence of this transition amplitude as compared with that coming 
from the spatial integral in the right hand side of (|5]), and take (e + e~\t\i) oc ipi(r = 0) 
as predicted by the nonrelativistic model of qq states. Note that in the nonrelativistic 
approach the D-w&ve qq state does not contribute to the di-lepton production. For this 
reason, the absence of splitting the 25*- and .D-states in the oscillator model is not very 
important from the point of view of the e + e~ mass spectrum. 

The applicability of the full fledged couple-channel formalism depends on the two 
space-time scales: the formation length Lf, associated with the i — > k transitions, 

L f = T- = ^2^2 ~ ^2 « 0.25/m - (24) 



kik ntf — ml 111 p' ~ m p o lGeV 



and the coherence length 



2E E m 2 n o 



Lc= k^ = M^Tq! W °- 75/m TGeV ' Jpfq- (25) 

associated with the transition 7* — > i. Strictly speaking, evaluation of diffraction matrix 
from equation (fLTD and of excitation amplitudes from (|T^) in terms of the color dipole cross 
section only holds if Lf > Rn and L c > R^. For the interesting to us region Q 2 < 1 GeV 2 
and for the p°(lS)- and p'(2S')-mesons the full fledged coupled-channel effects develop only 
at E ~ 5 — 8 GeV. However, at lower energy the only change is the decoupling of excitation 
of higher excitations from the photon and of off-diagonal diffractive transitions to and from 
higher excitations and we can stretch the formalism even down to E — 2 GeV. At this 
energy we have a single-channel problem with 7* — > p° excitation followed by diagonal 
p°N scattering. For the evaluation of the relevant diagonal elastic rescatterings we only 
need the p°N total cross section and the color dipole value of {p°\d-jp\p°) does still a good 
job for the almost energy independent Pomeron contribution. Simultaneously, (p°|o"|7*) 
enters only as the overall normalization and whether it is evaluated from eq. (|l^) or within 
different approach does not affect nuclear modifications of the e + e~ mass spectrum. 



4 Numerical results 
4.1 The input parameters 

We have performed the numerical calculations for the energies E = 2, 5, 10, 20 GeV at 
Q 2 = and 1 GeV 2 for the target nuclei 9 Be, 56 Fe and 207 Pb. For the parameterization of 
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the diffractive matrix and of the photon and vector meson wave functions see section 3. 
For the nuclear matter density in the light target nucleus 9 Be we use the oscillator shell 
model with the oscillator frequency adjusted to reproduce the experimental value of the 
root-mean-square radius of the charge distribution (r 2 )lg e = 2.51 fm f4lfl . For the target 
nucleus 56 Fe the parameterization of the nuclear density by a sum of Gaussians from ref. 
41| was used. For 207 Pb we use the Wood-Saxon parameterization of the nuclear density 



with parameters borrowed from ET]. 



4.2 The basis of vector meson states 

As stated above, in our numerical calculations we include 4 transverse resonance states. 
For the number of the off-diagonal rescatterings of the qq state we take v — 2. We checked 
that in our kinematic region the contribution from higher excitations and higher order off- 
diagonal rescatterings can safely be neglected. Furthermore, an approximation of the first 
two states and one off-diagonal rescattering of the qq state is sufficient for all the practical 
purposes. Our principal interest is in the interplay of nuclear effects and interference of 
the p°(lS) and p'(2S) and we did not include the numerically smaller contributions from 
the uj- and 0-mesons and their excitations. These long-lived resonances decay for the 
most part outside the target nucleus and the nuclear effects do not modify considerably 
their shapes. 



4.3 The presentation of results 

Due to the above mentioned form factor effect connected with the longitudinal momentum 
transfer and the suppression of the p'{2S) production by the node effect the amplitude 
(H) decreases strongly with increase of the mass of the e + e~ pair. In order to facilitate 
graphical presentation of the results, following Ref. |TU] we use the scaled amplitude 



T'(E, M) = M 2 T(E, M) /A . (26) 
Since in the absence of absorption effects fl5|) gives the amplitude oc A, in (|26"D we also 



introduced the factor I /A. Our principal numerical results for the mass spectrum are 
shown in figs. 1-6 in the form of 



|r , |2oc M« Ac 



A 2 dM 2 dpl 



p±=o 



We focus on the mass range M < 1.75 GeV and include as final states the p°(lS) and 
p'(2S) vector mesons. 

In order to see better the resonance behaviour of production amplitudes we also show 
in figs. 7-9 the Argand plots for X". 

Now we shall comment on the salient features of the Q 2 - , energy and nuclear target 
dependence of these mass spectra. 
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4.4 The nucleon target 

The reference e + e~ mass spectrum for the proton target and the photon energy E = 5 
GeV is shown by thick solid curves in figs. lb-6b. In the approximation of the energy 
independent dipole cross section a(p) it does not depend on the photon energy E. The 
mass spectrum exhibits the well separated p°(760) and p'(1480) resonance peaks. Recall 
the factor M 4 which makes the p'(2S) tail about flat at large M 2 . 

The suppression by the node effect is lifted with increasing Q 2 , see eq. (|21~1), and the 
comparison of the mass spectra for real, Q 2 = 0, and virtual, Q 2 = 1 GeV 2 , photopro- 
duction shows clearly the predicted rise of the p'(2S)/ p°(lS) ratio with increasing Q 2 
IE, 2~A | . This rise of the p'(2S) signal with Q 2 is clearly seen from a comparison of 



the Argand diagrams of fig. 7b and fig. 7b. 

4.5 Mass shift vs. nuclear form factor effects: the p° region at 
low energy 

In figs. la-6a we show the mass spectrum for low energy, E = 2 GeV, at which the p'{2S) 
production is negligible and the single-channel Glauber approximation (|T0D holds. The 
energy E = 2 is very low such that the p° would be expected to decay within nucleus. 
Although this reasoning is correct, our results show that the effect of the mass shift ([[]) 
due to the real part of the forward p°N scattering amplitude proves to be marginal in the 
nuclear modification of the inside component of the production amplitude. The point is 
that mass shift enters the exponent of the integrand in ([10|) via the extra phase 

= -t(b, z, zx)lva(j hh m - ■ r w m A {z - z x )J^(z - zi) (27) 

which must be compared to the phase fc 7 *^(z— Z\) = —{z — zi)/L c from the finite coherence 
length. Evidently, the significance of the mass shift for the inside component of the 
dilepton production amplitude is controlled by the parameter 



V = ~- r w ,n M f^L c » Q.l m ? 9 4/— g-. (28) 
' 2 P Is c M 2 + Q 2 \ 1 GeV V ; 

In the numerical evaluation in (|25J) we used the standard parameters of the p°N interaction 
as cited in section 3 and normal nuclear density. We see that the mass shift ([[]) amounts 
to a renormalization 

k rh -> k Th (l - rj) . (29) 

At E = 2 GeV the parameter r\ <C 1. On the other hand, the coherence length ( p5|) is 
very short, 

m 2 

L C (E = 2 GeV) = 1.5 fm • p 



M 2 + Q 2 ' 

much smaller than nuclear radii for heavy nuclei and even for the light Be nucleus it is 
comparable to the nuclear radius. For this reason nuclear effects are dominated by the 
attenuation of the p° due to the diagonal p°N p°N transitions and, most significantly, 
by distortions due to the nuclear form factor effects. 
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The results for the mass spectrum are shown in figs. la-6a by a solid curve. For the 
Be target the principal effect is the mass spectrum drops at M > 1 GeV much faster than 
for the free nucleon target. For the heavy Pb target the form factor oscillations lead to an 
effective splitting of the p° peak - the mass spectrum develops a minimum at M ~ m p o. 
At Q 2 = 1 GeV 2 the coherence length becomes still smaller and the distortion of the e + e~ 
mass spectrum by form factor effects become much stronger, such that the dip at the p° 
mass evolves already for the Be target. The corresponding Argand diagrams of figs. 8a 
and 9a span the mass range 0.5 < M < 1 GeV and exhibit a structure more complex 
than a single resonance loop. Our findings for real photoproduction off Fe and Pb targets 
are similar to those obtained in ref. for the target mass number A = 50 and A = 200 



in the approximation of uniform nuclear density. 

Notice that nowhere the splitting of the p° mass spectrum looks as a superposition 
of two Breit-Wigner peaks with the in-vacuum p° mass and the in-medium mass shifted 
by m 50 MeV as cited in section 3. The weak impact of the in-medium shift ([!]) on the 
p° splitting is obvious from dashed curves in figs. la-6a which show the mass spectrum 
obtained when the mass shift (lip is neglected, i.e., putting Re/Im = 0. The distortions 
of the mass spectrum change little, as a matter of fact the splitting of the p° peak is 
even enhanced somewhat in conformity to the rescaling (p9|) . Despite the slow rise of the 
parameter 77 with energy, in the range E =2-20 GeV of the interest for experiments at the 
Jefferson laboratory r\ -C 1. Furthermore, at higher energies the contribution from the 
in-medium decays decreases and the overall effect of the mass shift becomes still weaker. 

4.6 Higher energies: opening of the p'(2S) channel. 

The new feature of high energy photoproduction is opening of the p'(2S) channel and 
coupled-channel effects. For the comparison purposes, we show by the thick solid curve 
in figs. lb-6b the mass spectrum for the free nucleon target. The full coupled-channel 
results for nuclear targets comprise the attenuation, nuclear form factor effects, the effects 
of p°(lS) <-> p'{2S) transitions including the multistep transitions and the resonance mass 
shift due to Re/Im ^ 0. The relative importance of these effects depends on the mass 
region, energy and the target mass number. 

The importance of the off-diagonal transitions in the target nucleus can be judged from 
the comparison of the full coupled- channel results (thin solid curve) with those from the 
diagonal approximation (dotted curve) in which off-diagonal rescatterings are switched 
off, i.e., v = 0, and only the direct 7* — > p°(lS) and 7* — > p'{2S) transitions followed by 
elastic p°N and p'(2S)N rescatterings are included. One can see that in the p°-meson 
mass region 0.5 < M < 1 GeV the effect of the off-diagonal rescatterings for both the 
values of Q 2 turns out to be small. Indeed, the direct 7* — > p°(lS) transition is strong 
whereas the both transitions in the off-diagonal sequence 7* — ► p'(2S) — > p°(lS) are 
weak, see the description of diffraction matrix in section 3. Hereafter in our discussion 
of the off-diagonal effect we focus on the p'(2S) contribution because the contribution of 
intermediate states heavier than p'(2S)-meson is still smaller. 

Our results for real photoproduction, Q 2 = 0, in the p°(lS) mass region 56 Fe and 



207t 



Pb targets are close to that of ref. [ 10 1 obtained for real photons for the nucleus mass 



numbers A = 50 and A = 200 in the approximation of uniform nuclear density. 
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At Q 2 = 1 GeV 2 the coherence length is short also for E = 5 GeV and distortions 
of the shape of the e + e~ mass spectrum in the p°(15')-meson mass region by the nuclear 
form factor as described in subsection 4.5 persist for E = 5 GeV as well. The dominance 
of the form factor effects is obvious from the fact that these distortions change little from 
the full coupled-channel to diagonal case. 

At higher masses the coupled channel effects become more important. Here the direct 
7* — > p'(2S) transition is weak and there is a strong interference with the off-diagonal 
sequential transition 7* — > p°(lS) — > p'{2S) which contains a comparably weak transition 
p°(lS) — > p'{2S). As a matter of fact, for heavy target nuclei the sequential mechanism 
7* — > p°(lS) — > p'(2S) is found to dominate over the direct mechanism 7* — > p'(2S), 
which can be seen as follows. In figs. 1-6 we show by the dashed and long-dashed curve the 
pure p°(lS) and the p'(2S) contributions evaluated in the full coupled-channel approach, 
v = 2. The dot-dashed curve shows the pure p'{2S) contribution evaluated in the diagonal 
approximation, v = 0, and in this approximation the p'(2S) signal is much weaker than 
in the coupled-channel case in which the p'{2S) is fed by sequential transitions. This is a 
very interesting example when the p'{2S) production off heavy nuclei opens a possibility 
for extracting the matrix element (p'(2S)\a\p°(lS)) from the experimental data on the 
e + e~ mass spectrum in the p'(2S) mass region. The measurement of this matrix element 
could give a unique information on the overlap of the p°(lS) and p'{2S) wave functions. 

The p°(lS)-p'(2S) interference changes from destructive at masses M below the p'{2S) 
peak to constructive at, and above, the p'{2S) peak, reflecting the mass dependence of the 
relative phase of the p°(lS) — > e + e~ and p'(2S) — > e + e~ Breit-Wigner amplitudes: i) the 
mass spectrum develops a dip in between the p°(lS) and p'{2S) resonance peaks, where 
the coupled-channel, thin solid, curve goes below the thin solid curve for the pure p°(lS) 
contribution, ii) at, and beyond, the p'{2S) peak the coupled-channel, thin solid curve 
goes well above the long-dashed curve for the pure p'{2S) contribution, which testifies 
to a substantial contribution form the large-mass tail of the p°(lS) in the p'(2S) region. 
Evidently, that makes impossible an experimental extraction of the cross section of p'{2S) 
production in a probabilistic approach when the interference with the Breit-Wigner tail 
of p°(lS')-meson is neglected. 

4.7 The form factor effects in the p'(2S) region. 

For moderate energies and heavier nuclei, and also for larger Q 2 , we encounter the situation 
when the coherence length L c for the p'{2S) mass region becomes comparable to and/or 
smaller than the nuclear radius Ra- In this case the p'{2S) signal will be subject to 
distortions by the form factor effect in precisely the same manner as the p° signal at lower 
energies, see section 4.5. For instance, in real photoproduction off 56 Fe target at E = 5 
GeV the p'{2S) peak splits into two bumps with the dip at M m p /, see fig. 3b. The 
point that this dip is connected with the form factor effects is evident from the shift of 
the dip towards the smaller values of M, and the development of the secondary dip, with 
the increase of Q 2 , compare figs. 3b and 4b. The results for heavy target, 207 Pb, fig. 5b 
to fig. 5c to fig 5d show clearly how the dip-bump structure moves to higher masses M 
with the increase of the coherence length L c with rising energy E. 

A comparison of Argand diagrams for the free nucleon target, fig. 7, and nuclear tar- 
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gets, figs. 8 and 9, shows that the resonance loop corresponding to the p'(2S')-meson 
becomes well visible only at E :> 10 GeV. However, even in this energy region the contri- 
bution to the amplitude from the Breit-Wigner tail of p°(lS) meson cannot be neglected 
as compared with that from the p'(2S) meson. 



5 Conclusions 

We have performed a coupled-channel analysis of nuclear-medium modification of the 
p°(lS)- and p'(2S')-mesons in the coherent 7* — > p°(lS)A, p'(2S)A — > e + e~A reaction in 
the kinematic region E ~ 2 — 20 GeV and Q 2 < 1 GeV 2 . Our findings on the interplay of 
the inside and outside decays, in-medium modifications of the inside component and the 
coupled-channel effects can be summarized as follows: 

(i) In the p°(lS) meson mass region 0.5 < M < 1 GeV the effect of the off-diagonal 
rescatterings is small. For heavy nuclei in this mass region our results agree with those 



obtained in ref. |L0| within the one-channel approximation. 

(ii) The off-diagonal rescatterings become important at M > 1 GeV. The main off- 
diagonal contribution is the p'{2S) production through the sequential mechanism 7* — > 
P°(1S) — > p'(2S). This mechanism dominates the cross section of p'(2S) production for 
heavy nuclei. 

(iii) At low energies E < 5 GeV the shapes of the p°(lS) and p'{2S) resonances are 
strongly affected by the nuclear effects connected with the interference interplay of the 
resonance decays inside and outside the target nucleus and form factor effect connected 
with the longitudinal momentum transfer. 

(iv) The p'{2S) resonance is seen well only at E > 10 GeV. Even at high energy its shape 
is affected by the interference with the Breit-Wigner tail of p°(15')-meson which must be 
included properly in an analysis of the experimental data on the e + e~ mass spectrum in 
the p'{2S) mass region. 



In our analysis we focused on the dilepton decay mode. Evidently, very similar effects 
must be observed in measurements of the p°(lS) and p'{2S) production through the im 
decay mode. The major difference from the dilepton mode would come from to the final- 
state interaction of the irn system, which would reduce the relative contribution of the 
inside component as compared with that for the e + e~ mode. For this reason a comparison 
of the mass spectra for these two cases would be of great interest. It is especially interesting 
at low energies E ~ 2 — 5 GeV for heavy nuclei when the inside component is large for 
the e + e~ decay mode and will be strongly suppressed by the final-state absorption for the 
7T7T mode. The comparative theoretical analysis of the coherent p°(lS) and p'{2S) photo- 
and electroproduction for the e + e~ and mi decay modes is now in progress. 
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Figure captions: 



Fig. 1 - The e + e~,p + p~ rescaled mass spectrum (in arbitrary units) for coherent real 
photoproduction, Q 2 = 0, off the 9 Be nucleus for the incident beam energy E = 
2, 5, 10 and 20 GeV. The legend of curves in box (a): The solid curve is a prediction 
of the Glauber approximation with mass shift flTU| ) from the reggeon amplitude Q2"3| ) 
included, the dashed line shows the results obtained neglecting the mass shift. The 
legend of curves in box (b): The thick solid curve is the spectrum for the nucleon 
target, all other curves are for nuclear target. The thin solid curve is the full coupled- 
channel calculation, v = 2 in expansion (Ip. The dotted curve is a result from the 



diagonal approximation, v = in expansion flltl), the p°(lS)-p'(2S) interference in 
the e + e~ and/or p + p~ decay channels included. The dashed curve shows the pure 
p° signal in the coupled-channel calculation, v = 2. The long-dashed curve is the 
p'{2S) signal in the coupled- channel calculation, v = 2. The dot-dashed curve is 
the p'(2S) signal in the diagonal approximation, v = 2. Boxes (c) and (d) are the 
same as box (b) but for E = 10 and 20 GeV, respectively, and with omission of 
mass spectrum for the free nucleon target. 

Fig. 2 - The same as Fig. 1 but for virtual photoproduction, Q 2 = 1 GeV 2 . 

Fig. 3 - The same as Fig. 1 but for the target nucleus 56 Fe. 

Fig. 4 - The same as Fig. 2 but for the target nucleus 56 Fe. 

Fig. 5 - The same as Fig. 1 but for the target nucleus 207 Pb. 

Fig. 6 - The same as Fig. 2 but for the target nucleus 207 Pb. 

Fig. 7 - The Argand plots for the scaled amplitude V = M 2 T(E, M) for (a) real, Q 2 = 0, 
and (b) virtual, Q 2 = 1 GeV 2 , photoproduction off the free nucleon target. 

Fig. 8 - - The Argand plots for the scaled amplitude ( p6|) at Q 2 = 0. The solid, dotted, 
and dashed curves show the results for the target nucleus 207 Pb, 56 Fe, and 9 Be, 
respectively. The curves are given in arbitrary units. The mass intervals are 0.5-1 
GeV for E = 2 GeV, and 0.5-2 GeV for E = 5, 10, and 20 GeV. The spacing of mass 
points along the curves is 0.25 GeV, the arrows show the direction of the increasing 
mass. 

Fig. 9 - The same as Fig. 8 but for Q 2 = 1 GeV 2 . 
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